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ABSTRACT 



12 C/ 13 C isotopologue abundance anomalies have long been predicted for gas-phase 
chemistry in molecules other than CO and have recently been observed in the Taurus 
molecular cloud in several species hosting more than one carbon atom, i.e. CCH, CCS, 
CCCS and HC 3 N. Here we work to ascertain whether these isotopologic anomalies 
actually result from the predicted depletion of the 13 C + ion in an oxygen-rich optically- 
shielded dense gas, or from some other more particular mechanism or mechanisms. 
We observed A3mm emission from carbon, sulfur and nitrogen-bearing isotopologues 
of HNC, CS and H 2 CS at three positions in Taurus (TMC1, L1527 and the NH 3 peak) 
using the ARO 12m telescope. We saw no evidence of 12 C/ 13 C anomalies in our ob- 
servations. Although the pool of C + is likely to be depleted in 13 C, 13 C is not depleted 
in the general pool of carbon outside CO, which probably exists mostly in the form of 
C°. The observed isotopologic abundance anomalies are peculiar to those species in 
which they are found. 

Subject headings: astrochemistry . ISM: molecules . ISM: clouds. Galaxy 
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1. Introduction 



Observations of isotopically-substituted versions of interstellar molecules - their so-called 
isotopologues - have become the usual way of deducing atomic isotopic abundance ratios for 
common products of stell ar nucleosynthe sis - C, N, O, S, Si , CI etc in the interstellar medium 



(ISM), as summarized by 



Wannied (11980) and 



Wilson! (| 19990 . However, the relative abundances 



of particular isotopo logues can differ f r om those inhe r ent in the gas at lar ge owing to peculiarities 



of their chemistry (|Watson et al. 



1976 



Watson 



1977 



Langer et al 



1984), and recognizing or 



explaining anomalous isotopologic abundance ratios is an ongoing challenge for astrochemistry. 

The most commonly-observed example of this phenomenon is with CO in diffuse clouds, 
where the the CO/ 13 CO ratio may be larger or smaller than the inherent C/ 13 C ratio in the 



ambient gas depending on whether selecti ve photodiss ociation (Bal 



chemical fractionation (Watson et al. 



S onnentrucker et al. 



2007 



Burgh et al. 



1976)) dominates (Wilson et al 



2007 



Sheffer et al 



2008 



v & L anger 



Liszt 



992; 



2007; 



1982)or 



Liszt & Lucas 



Visser et al 



1998; 



2009) 



Selective photodissociation (increasing CO/ 13 CO) occurs in gas seen toward a few hot stars while 
chemical fractionation (behaving oppositely) is a more general phenomenon that is enhanced at 
somewhat highe r CO column densitie s. Chemical fractionation is also now seen in the envelopes 



of dark clouds (Goldsmith et al. 



2008). 



Gas-phase molecular cloud chemistry actually makes a rather general prediction about 
isotopologic abundances in strongly shielded regions: CO and molecules that form directly from 
it (HCO + , H 2 CO and CH 3 OH) should show a common C/ 13 C ratio that reflects the composition 
of the gas reservoir, while other sp ecies that form from the pool of free carbon outside CO 



should be strongly depleted in 13 C (|Watson 



1977 



Liszt 



1978 



Langer et al 



1984). This occurs 



'Based on observations obtained with the ARO Kitt Peak 12m telescope. 



because CO is the repository of the overwhelming majority of the gas-phase carboro and carbon 
is liberated from CO mainly as the dissociation product of destruction of CO by a small quantity 
of cosmic-ray ionized He + (i.e. He + + CO — » C + + O + He). The resulting C + ions quickly 
interact with the ambient CO, and the fact that 13 CO is more strongly bound than 12 CO results in a 
pref erential deposition of 13 C into 13 CO through the reaction 13 C + + 12 CO -> 12 C + + 13 CO + 34.8 



K (Watson et al. 



(119761) . see Appendix A here). In this way 13 C + preferentially disappears from 
the pool of C + and from the general pool of carbon available to form most species. The C + / 13 C + 
ratio increases by a factor approaching exp(34.8/T K ) but the CO/ 13 CO ratio changes little because 
the C + /CO abundance ratio is so small, below 10~ 3 . 

It quickly became apparent that there was some tension between this prediction and 
observations of interstellar gas showing C/ 13 C ratios that were, if anything, below the terrestrial 
value of C/ 13 C = 89. Now it is ack nowledged tha t the local interstellar carbon isotope ratio is 



in the range 60-70, but at the time, IWatsonl (|1977h suggested a mechanism involving selective 
depletion of carbon onto grains, offering the possibility of lowering the C/ 13 C ratio. Detailed 
time-dependent models showed v ery complex behaviour in the isotopologic ratios but largely 



failed to bear out this suggestion (ILiszt 



19781) . finding that 13 C enhancement outside CO only 



occured as molecules disappeared from the gas. 

Because of the possibility of fractionation in fully molecular gas, some observers changed 
the focus of their work to concentrate on measurement of carbo n isotope ratios using only 



19931) but it seems fair to say that 



the isotopologues of carbon monoxide (ILanger & Penziasl 

anomalously large C/ 13 C ratios have never been se en in the molecular cloud- H II region complexes 



that are typically used for isotope determinations (IWanniei 



1980; 



Wilson 



1 



1999). Indee d, a recent 



surve y of the carbon iso tope ratio in CN does not seem to have been affected (IMilam et al. 



and 



Tercero et al. 



2005) 



(1201 Oh report C/ 13 C = 45+ 20, S/ 34 S= 20 + 6 from a survey of sulfur-bearing 



/34 



2 at least as long as [O] > [C] 
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carbon chain molecules in Ori KL. Fractionation effects may be minimized at the somewhat 
higher temperatures of clouds near H II regions but are not entirely eliminated. The absence of 
observable effects from the predicted carbon fractionation chemistry has been a lingering, if not 
entirely obvious, mystery. 

By contrast, anomalies of two kinds have now been noted in cold gas in/near the 



cyanopoly yne peak in t 



le Taurus Molecular C 



atoms, and 



Sakai et al. 



(120071) and 



Sakai et al. 



oud ( TMC) for species contai ning several carbon 



(|2010 ) appealed to the work of 



Langer et al 



(1984) 



for an explanation. C is strongly 



isotopologues of CCH (|Sakai et al. 



put ve ry unequally lacking in both of the singly-substituted C 



20101) and is deficient in 13 CCS but not C 13 CS. A small effect 



may also be present in HC 3 N in Taurus (iTakano et al. 



1998) and 



Martin et al. 



(12010 ) recently 



reported N(CCH)/N( 13 CCH) > 1 10 (3<r) in M 82. 

The gas-phase carbon isotopic chemistry is the subject of this work. Motivated by the 
apparent contradictions among the model predictions for a strong depl etion of 13 C, the ab sence of 



Milam et al. 



(2005) for CN 



any obvious effect in subsequent measurements like those reported by 
and 13 CN, and the rec ent invocation of a gene ral 13 C depletion to explain anomalies seen in the 



nearby Taurus Cloud (ISakai et al. 



2007, 



2010(), we observed isotopologues of species containing 



a single carbon (CS, HNC and H 2 CS) toward three positions in Taurus, as described in Sect. 2. 
Section 3 summarizes the observational situation including our new measurements. Having found 
no evidence that 13 C is depleted in the species we observed, and probably in the general pool of 
carbon in the gas, we discuss the overall carbon chemistry in Sect. 4, inquiring how 13 C depletion 
can be avoided and the implications for the gas phase chemistry if it does not actually occur. 
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2. Observations and conventions 

The new observations reported here were taken at the ARO 12m telescope in 201 1 March and 
April. We observed HNC, CS, H 2 CS and their singly-substituted 13 C, 15 N and 34 S isotopologues 
at three positions in Taurus: the TMC1 cyanopolyyne peak (a = 4 h 41 m 42.88 s , <5=25°41'27"), 
L1527 (a = 4 h 39 m 53.89 s , 5=26°03'11") and the NH 3 peak (a = 4 h 41 m 23 s , £=25°48'13.3"; all 
coordinates are J2000). All the data were taken by position-switching to an off-position 8.5' NE 
of the cyanopolyyne peak, using the dual-polarization mode of the autocorrelator configured for 
24.4 kHz channel spacing (0.073 kms 1 at 100 GHz) and 48.8 kHz spectral resolution. Typical 
system temperatures were 180 - 240 K depending on line frequency and source elevation. All 
velocities are given with respect to the kinematic local standard of rest and the line temperatures 
are given in terms of the native scale at the ARO 12m, T^« 0.85 T mb . Line frequencies were taken 
from the online splatalog (http://www.splatalogue.net) and are given in Table 1 representing the 
new results. 

In this work we do not call out the most-abundant isotope explicitly: C implies 12 C. The 
integrated intensity of spectral lines for species X, in units of K km s _1 , is denoted by W(X). 



3. Observational results 
3.1. New results 

Spectra are shown in Fig. 1 and numerical results are summarized in Table 1 giving profile 
integrals and Table 2 giving ratios of line profile integrals and implied carbon isotope ratios, 
assuming values for N/ 15 N or S/ 34 S as indicated. 

The main isotopologues of HNC and CS are optically thick and self-absorbed by weakly- 
excited molecules within the foreground envelopes; it is not possible to measure the C/ 13 C ratios 
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Fig. 1. — Observed line profiles toward three positions in TMC1 
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in these species directly and we must assume values for the isotopologic ratios in sulfur and 
nitrogen. For nitrogen the interstellar N/ 15 N ratio is now seen to be only slightly larger than 



the terrestrial ratio of 270. 



Lisetal 



(1201 Oh found N/ 15 N = 334+50 (3cr) toward Barnard 1 and 
Adande and Ziurys (in press) find N/ 15 N= 290 +40 at the Solar Circ le from a large-scale galactic 



survey of CN emission. In local diffuse clouds 



Lucas & Lisztl (11998b found N/ 15 N=240+25. The 



range of values implied for HNC/HN 13 C in Table 2 corresponds to HNC/H 15 NC = 250-330 and 
the statistical uncertainties are quoted at each end of the range. 



No systematic deviations from the S olar System va 



1980; 



ue of S/ 34 S= 22.7 have been see n in 



Wilson 



1999; 



Tercero et al. 



2010) and 



specific isotope studies of the local ISM (Wannierj . 
the implied values of CS/ 13 CS quoted in Table 2 correspond to this value and the statistical 
uncertainties of the current measurements. Our measurements of H 2 CS isotopologues toward 
TMC1 are in accord with this. Some uncertainty accrues from the unmeasured optical depth of 
the rarer isotopologic lines as discussed in the following subsection. 

The result for the isotopologic carbon ratio in H 2 CS is of marginal significance and will not 
be discussed further but it seems clear from Table 2 that there is no strong depletion of 13 C in 
HNC or CS. 



3.2. Correction for finite optical depth 



3.2.1. CS 



The presence of a finite optical depth in the C S line would cause compression of the 



W(C S)/W( CS) intensity ratio, resulting in an und erestimate o f the 



N(CS)/N( 13 CS) column density ratios. For TMC1, 



Pratap et al 



N(C 34 S)/N( 13 CS) and 



(1997) derived N(C 34 S) 



2.0 x 10 12 cm" 2 , consistent with an estimate N(C 34 S) = N(CS)/22.7 = 5 x 10 13 cm" 2 /22.7 



2.2 x 10 cm from the work of Oh ishi et al. 



(1992). ForN(C 34 S) = 2.0 x 10 12 cm" 2 and peak line 
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brightness of 0.72 K - 0.85 K, depending on beam efficiency, we calculate that the central optical 
depth of the C 34 S J=2- 1 line toward TMC1 must be below 0.15 and the excitation temperature 
in the lower part of the rotation ladder must be fairly high, above about 8 K for J=2-l, based on 
modelling of the excitation. In this case the 13 CS line is quite optically thin and the observed 
isotopologic intensity ratio is at most 5% below the intrinsic abundance ratio. Accordingly, the 
entry corrected for optical depth in Table 2 is the uncorrected value increased by 5%. 

Estimates of the optical dept h correction tow ard the other pointing positions do not seem 
possible. Toward the NH 3 peak LPratap et al.l (|1997|) derived a higher column density N(C 34 S) = 
4.3 x 10 12 cm -2 while the observed line in this work is only slightly broader and noticeably less 
bright. Excitation solutions for such high column density in the face of relatively weak emission 
indicate high optical depth at lower excitation temperatures (5 K) that do not seem consistent with 



the physica 



Pratap et al 



condit ions toward the NH 3 peak, that are not less dense or colder than toward TMC1. 



(|1997|) cited problems with their modelling of the CS lines, such that they were 
compelled to use a very small isotopologic abundance ratio CS/C 34 S = 12-14. Although they cited 
a reference to the isotopic composition of cosmic rays in defense of this assumption, it is without 
precedent in molecular gas and would lead to improbably small values of the CS/ 13 CS ratio if used 
in this work. We conclude that we cannot make a reliable correction for optical depth correction 
except toward TMC1 where all indications are that the opacities in the rarer CS isotopologues are 
small. 



3.2.2. HNC 

For HNC, the 13 C-bearing isotopologue is the more abundant of the two variants used to 
derive the C/ 13 C ratio, so the effect of optical depth is opposite to that in the case of CS: the 
derived HNC/HN 13 C ratio decreases with increasing optical depth in the rare isotopologues. 
Indeed, the optical depths encountered in HNC are higher than for CS and the correction for 
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finite optical depth is rather larger. Unfortunately we did not resolve the hyp erfine structu r e of 



HN 13 Cand cannot duplicate the excitation temperature/opacity derivation of 



Padovani et al. 



(1201 lh 



who observed at other positions in Taurus having higher optical depth than occurs toward TMC1. 



Pratap et al. 



(119971) gave N(HN 13 C) = 4.2 x 10 12 cm" 2 toward TMC1 while a slightly smaller 



uppe r limit N(HN 13 C) < 2 x 10 14 cm 2 /60 < 3.3 x 10 12 cm 2 can be inferred from the tabulation 



of 



Ohishi et al 



(|1992|) . For the comparatively weak emission lines we observed, these column 
densities imply J= 1-0 excitation temperatures around 4 K and optical depths of order 2 in HN 13 C, 
so that even such a rare isotopolgue would be somewhat self-absorbed. With optical depths 
this large, simply requiring that HNC/HN 13 C be above 60 limits the possible range of assumed 
values for N(HNC)/N(H 15 NC) to be nearer the large end of the likely range. However, the small 
excitation temperatures implied by our observed brightnesses and the higher of the previously 
measured column densities are not consistent with the densities that are inferred for the Taurus 
cores, which imply higher excitation temperatures, lower optical depths and smaller column 
densities. 



We conclude that it is not possible to derive an entirely accurate carbon isotopologue ratio 
in HNC based on correcting the observed intensity ratios for optical depth but the errors do go 
in a known direction. Table 2 shows a range of carbon isotope ratios derived from HNC for an 
optical depth 0.6 in the HN 13 C line toward TMC1, correspon ding to an excitatio n solution with 



N(HN 13 C) = 2 x 10 12 cm 2 that is consistent with the work of 
understimate the optical depths and so overestimate the C/ 13 C ratio. 



Ohishi et al. 



(1992) . This could still 



3.3. Previously-reported 13 C anomalies 



These are summarized in Table 3 where all the results are cast in terms of actual or implied 
C/ 13 C ratios, using the terrestrial sulfur isotopic ratio as appropriate. There are two effects 
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that must be considered: 13 C is strongly lacking in some isotopologues but to different degrees 
depending upon placement. Thus both 13 CCH and C 13 CH are underabundant but moreso in the 
former (a factor 1 .6) and in CCS only the more tightly-bound isotopolgue 13 CCS exhibits the 
effect; in each case, the isotopolgue with an inboard 13 C is more abundant. In C3S only 13 CCCS 
was observed, but is heavily deficient. In terms of chemical models these two considerations are 
separately reflected in, on the one hand, the overall C/ 13 C ratio in the gas, and, on the other, the 
detailed path to formation (and perhaps in situ fractionation after formation) of the individual 
species studied. 



Sakai et al. 



(120071) and 



Sakai et al. 



(120 lCy discussed in exquisite detail the formation routes to 
CCS and CCH, concluding that the unequal depletions of the 13 C-bearing isotopologues were the 
result of formation reactions in which the carbon atoms were not equivalent (for CCH, C + CH 2 
— > CCH + H and not C 2 H 3 + + e — » CCH + H 2 ), rather than alterations that might have occurred 
after formation (for instance 13 CCH + H -> C 13 CH + H + 8.1 K or 13 CCS + S -> C 13 CS + S + 
15 K). They also co ncluded that they had confirmed observationally the phenomena modelled by 



Langer et al. 



(|1984f) whereby the relatively famous fractionation reaction involving CO and C + 
induces a very general vanishing of 13 C from the gas at large as discussed in the Introduction here. 



4. Carbon isotopic and isotopologic ratios in the Taurus cores 



The C/ 13 C ratios inferred from our work are consistent with other measurements of the 
carbon isotopic ratio in the nearby ISM that find C/ 13 C = 60-70 and show no sign of the 13 C 
depletion that is predicted for the gas phase chemistry. The CS isotopologues we observed would 
not serve as the basi s for substantiall y depleted CCS but are instead close to the CCS/C 13 CS ratio 



54 + 2 discussed by 



Sakai et al. 



(|2007l) . suggesting that CCS could indeed form directly through a 



reaction like CH + CS — » CCS + H as hypothesized by lSakai et al. 



(2007). 
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In fact, Furuya et al.l (|201lf) recently reached something of the same conclusion through 



chemical modelling, without the new observational result for CS. They consi dered a time 



dependent model for TMC1 with all the carbon initially in the form of C + (see 



Lisztl (12009)), in 



which the abundance of CO is relatively small even when it begins to deplete out of the gas onto 

grains; the CO column in their model is al so rather small compared with N(CO) = 500 N(C 18 0) 

I J 

= 1.8 x 10 18 cm" 2 from lPratap et al.l (|1997|) . Because of the assumed initial conditions a very large 
proportion of the gas-phase carbon remains outside CO, and, for reasons that are discussed in Sect. 
5, 13 C depletion in the gas as a whole is generally modest. C/ 13 C is typically of order 1 10 after the 
earliest times. In this case the hypothesis of a rearrangement effected by reaction of CCS and H 
yielded CCS/ 13 CCS = 230, CCS/C 13 CS = 54, as observed. That is, the rare isotopologues both 
form with C/ 13 C = 1 10 but the hypothesized reaction C 13 CS + H -> 13 CCS + 17.4 K is strongly 
energetically favored at 10 K driving CCS/ 13 CCS up while CCS/C 13 CS falls to half the C/ 13 C ratio 
in the gas (110/2). Our results suggest that 13 CCS and C 13 CS actually form with CCS/ 13 CCS and 
CCS/C 13 CS =* 70. 



Furuya et al.l (|2011r) also showed that the known rearrangement reaction for CCH with 
ambient H (that has a smaller exothermicity 8 K) could yield unequal isotopologic abundances to 
the degree observed (a factor 250/170) but the relatively small amount of 13 C depletion in their 
model did not come close to reproducing the very large individual CCH/ 13 CCH or CCH/C 13 CH 
ratios that are observed (170 and 250; Table 3). Even so, their model generally predicts larger 
values for N(CS)/N( 13 CS) than we observe. It is not obvious that this scheme would function 
as well in a gas that was even more weakly fractionated, if it caused very small values for the 
CCS/ 13 CCS rato. 

The smaller isotopologic discrepancies that are observed in HC3N can likely be achieved 
by relatively weakly-endothermic rearrangement reactions with hydrogen in an unfractionated 
gas, consistent with the lack of fractionation we infer for HNC. In summary, the observational 
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situation is complex and somewhat mysterious. CS, HNC and HC3N suggest that there is little 
fractionation in the carbon pool at large while the seeming absence of fractionation in CCS/C 13 CS 
may be accidental and the anomalies in CCH remain unexplained. Disparities among the carbon 
isotopologues in individual species with more than one carbon atom seem best explained by 
rearrangement reactions with atomic hydrogen but the very large isotopologic ratios seen in CCH 
remain beyond the reach of current explanations and the disparity that is observed in CCS may be 
inconsistent with rearrangement in a gas that is too-lightly fractionated. 



5. Where is the CR ionization-driven fractionation? 

The possibility of a strong overall carbon fractionation in species other than CO has been 
recognized ever since the existence of the C + + CO fractionation reaction was first remarked. 
Nonetheless, it was not observed in surveys that sought to determine the carbon isotope ratio but 
found the same values in CO and CS and CN etc, even though the latter species do not form from 
CO and could have shown strong fractionation effects even in relatively warm gas. Even now, 
despite the carbon isotopologic anomalies that are seen at very low temperatures (10 K) in Taurus 
in some polyatomics bearing more than one carbon, the gas there can not obviously be said to be 
fractionated as a whole. Here we ask why this is so. 



5.1. The role of He 



The prediction of strong overall fractionation arises in models of the gas-phase chemistry of 
dense dark ox ygen-rich material (O/C > 1) where CO is by far the largest repository of gas-phase 



carbon atoms (ILanger et al 



1984|) . A few molecular species form directly from CO - HCO + in the 
gas, H 2 CO and CH 3 OH from successive hydrogenation of CO on grains - and will share its C/ 13 C 



ratio. Other carbon-bearing trace molecules in the gas result from the relatively small amount of 
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carbon that is liberated when He + , ionized by cosmic rays, slices CO apart into C + and O via the 
reaction He + + CO — » C + + O + He; the reaction rate constant is 1.6 x 10 9 cm 3 s _1 (see Table 
4 for important reactions and rate constants). The only reaction competing effectively with this 
one for the attentions of He + is the charge-exchange ionization of H 2 , He + + H 2 — » He + H 2 + , 



whose currently-accepted reaction rate constant is 7.2 x 10 15 cm 3 s 1 (Table 4). 



Langer et al 



(| 19841) included instead the reaction He + + H 2 — > H + + H + He with the much larger rate constant 
1.5x10" 13 cm 3 s _1 , which significantly dampened the effects of fractionation compared to what 
would be derived in more recent chemical models: that reaction is now considered to be negligibly 
slow in cold gas, see Table 4. 

The ratio of reaction rates and abundances determines whether CO or H 2 controls the 
abundance of He + but for X(CO) > 4.5 x 10" 6 , when CO dominates using currently-accepted 
reaction rate constants, every cosmic -ray ionization of an He atom results in the production of a C + 
ion. Because He is so much more abundant than carbon (i.e. CO), the volume rate of production 
of C + ions and the rate at which carbon is liberated from CO is approximately 1000 times higher 
than the rate at which cosmic rays interact with CO directly; this factor is independent of the 
cosmic ray ionization rate n(He + ) and independent of X(CO) when the reaction with CO is the 
main destroyer of He + . 

It is important to stress just how rapidly the chemical ionization of carbon occurs via 
He + + CO, because the only competing mechanism for liberating atomic carbon from CO 
is photodissociation. In quiescent, shielded, cosmic-ray heated gas, cosmic-ray induced 
photodissociation typically occurs at rates that are a much smaller multiple of the cosmic ray 
ionization rate, and is some 50 times slower than He + + CO. The point is that in these models it is 
difficult to contrive to liberate carbon from CO except, initially, as C + , leading to fractionation. 



5.2. The role of the CO + C + fractionation reaction 



Just as He + reacts most commonly with CO, so does C + , although in this case there is a 
wider range of competitive reactions (see Table 4). But to the extent that C + interacts with CO, 
it preferentially deposits 13 C + back in CO. Meanwhile, the abundance of C + is so small relative 
to that of CO (typically less than 0.1%) that any effect on the isotopologic abundance ratio in 
CO is imperceptible. The final result is that when free atomic carbon is primarily liberated and 
maintained in the form of C + , the pool of gas-phase carbon outside CO inevitably becomes 
depleted in 13 C + , hence in 13 C, as long as CO remains sufficiently abundant in the gas to dominate 
the neutralization of He + . CO and the free carbon pool coexist but with two distinct C/ 13 C ratios 
that are separately passed on to their descendant molecules. 



5.3. So indeed, why is a more general 13 C depletion not seen in Taurus? 



To blunt the effect of C + fractionation it is necessary to interfere with the C + -CO interaction 
and there are various choices for doing so, for instance; i) a low CO abunda nce, such as might 



Furuva et al. 



(2011) 



occur at later times after depletion onto dust (see also the early-time model of 
in which CO never attains its full abundance); ii) very large abundances of species like OH and 
O2 that react rapidly with C + ; iii) taking the high-metals case where sulfur and silicon are not 
assumed to be strongly depleted from the gas so th at they dominate the neutralization of C + 



through charge exchange; iv) the presence of PAH (|Wakelam & Herbst 



2008 



Lepp et al. 



1988) 



to neutralize C + more rapidly, although these are not always assumed to survive in dense gas. 
Nonetheless, fractionation at the level of a factor two or so persists in the C + until CO and C + are 
all but gone from the gas. 

Preventing fractionation of the pool of molecules outside the small family that forms from 
CO really requires that the reservoir of carbon outside CO must reside mostly in the form of 
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neutral atomic carbon with perhaps some small admixture of C + . Given that C° typically reacts 
with neutral molecules at least 20 times more slowly than does C + , the proportion of C° to C + 
should be very large and the real question is how this might come about. One possibility that 
is capable of maintaining very large free neutral atomic carbon fractions is the so-called high 
ionization phase of bistable solutions of the chemical equilibrium equations, whose relevance has 



recently been disc ussed in terms remarkably similar to the considerations here for contrq 



the C + abundance (Char nlev & Markwick 



2003 



Wakelam et al 



2006; 



Boger & Sternberg 



ling 



2006). 



Extension of the bistability discussion to include fractionation effects is clearly of great interest. 

It is generally not possible to create a sufficient pool of neutral carbon using just the in situ 
flux of CO-dissociating photons generated by t he cosmic ray-induced electronic excitation of H 2 



(Prasad & Tarafdar 



1983 : 



Sternberg et al. 



1987|) because these p hotons dissociate CO at rates that 



Langer et al. 



1984) assumed a 



are a relatively small multiple of the cosmic ray ionization rate: 
multiple of 10. As such, cosmic -ray induced photons liberate neutral atomic carbon much more 
slowly than CO is sliced into C + and O by He + . Because the naturally-arising local flux of CO 
photodissociating photons inside dark gas is so weak, some attempts to explain the existence of 
even a relatively small amou nt of C° in dense ga s invoke special mechanisms such as inversion of 



the normal [0]/[C] > 1 ratio (Langer et al. 



1984) . penetra tion of ambient uv 



heavily clumped medium and/or recent cloud formation ([Phillips & Huggins 



ight i nto a porous, 



19811), turbulent 



diffusion that cyclically exposes heavi 



cloud surface (Boland & de Jon 



1982; 



y-shielded gas to the ambient uv rad iation field near the 



Will acy et al 



2002; 



Xie et al 



the so-called high-ionization phase of bistable chemical reaction schemes (|Flower et al 



1995 ) and the presence of 



1994) 



The inference of a large pool of free atomic carbon in dense dark gas is really not an 
observational problem p er se because substa ntial columns of C° are typically found toward 



and ar ound dark clouds (IFrerking et al 



0.1 (Maezawa et al. 



1989h including toward TMC1 where N(C°)/N(CO) 



1999). This is very large compared to the values X(C + )/X(CO) 
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~ 3 x 10" 8 /8 x 10" 5 ~ 4 x 10" 4 that arise from the default parameters in our toy model (Appendix 
B). However, the observational situation is complicated by superposition of different density 
regimes along the line of sight and it is hard to assess how much of the observed atomic carbon 
actually exists inside the darkest gas. This is especially true given the recent revelation of the 



preponderance o f low-A v material in Taurus and other dark cloud com plexes (|Pineda et al 



Cambresv 



1999). Extrapolation of the results of Bensch et a 



suggests that perhaps 1-2% of the carbon budget (|Gerin et al. 



2010 : 



(2003) to much higher N(CO) 



2003b could exist as C° at N(CO) 



10 18 cm 2 . This would be about 50 times more carbon than exists in C + . 



5.4. Relevance to other environments 



The possibility of chemical fractionation in cold cores raises important questions for 
understanding carbon isotopic abundances in planetary and proto-planetary systems including the 
Solar System and proto-Solar nebula. For example, is it correct to ascribe the difference between 
the Solar ratio C/ 13 C = 89 and that in the nearby ISM (60-70) entirely t o chemical enrich ment 



since the birth of the Solar System or might fractionation be responsible (Smit h et al. 



20 lib ? Why 



are there differences between the carbon isotop e ratios measured in various comets and betwee n 



those measured in comets and that in the Sun (Crovisier et al 



2009 



Mumma & Charnley 



201 lh . 



Other questions arise in matters involving other elements, for instance the difference between 
the N/ 15 N ratio in the Earth (270) and Sun (420) and the wide disparities in D/H measurements 
between that in the earth's oceans (1.6 x 10~ 4 ), the more nearly cosmological ratios seen in the 
outer planets (3 -5 x 10~ 5 ) and the much higher than telluric D/H ratios seen in mo st comets 



(Crovisier et al 



201 lh . 



2009 



Mumma & Charnley 



201 lh . although, apparently, not all (IHartogh et al 



Many of the effects discusse d here occur prominently in recent models of fractionation 
chemistry in protoplanetary disks (IWoods & Willacyl 120091) . which illustrate the complexity of 
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relating isotope ratios in the proto-planetary nebula or planetary disk (when formed) to that in 
the ambient natal material. Fractionation varies with disk radius and disk height and evolves 
with time under the com bined influence of the proto-stellar and interstellar radiation fields. 



Woods & Willacy! (|2009t) conclude that Solar System cometary material has been reprocessed, 
raising the question of whether any memory of conditions in natal molecular material persists into 
fully-formed planetary systems. 



6. Summary 

Since the initial recognition of the carbon isotope fractionation reaction a conflict has existed 
between the very general prediction of a strong 13 C depletion in molecules other than CO and the 
general absence of observable effects in surveys of the C/ 13 C isotopic abundance ratio deduced 
form such common species as CS or CN. 13 C depletion is predicted to occur when carbon is 
liberated from CO by the reaction He + + CO — > C + + O + He and remains in the gas as C + to 
interact and fractionate with CO. Eventually, CO depletion onto grains will blunt the effect at late 
times but in the meanwhile quite large variations in the abundances of 13 C-bearing molecules are 
predicted for all species that do not form directly from CO. The point is that CO remains very 
nearly unfractionated as long as it is the main carbon reservoir in the gas and X(C + )/X(CO) is very 
small, so that species like HCO + , H 2 CO and CH 3 OH will also be unfractionated. But the relatively 
small pool of carbon that exists outside of CO to form other molecules than those in the CO family 
becomes strongly depleted in 13 C + and molecules that form from it (most carbon-bearing species) 
will also show strong 13 C depletion. 

However, strong anomalies in CCH and in some isotopologues of CCS and CCCS have 
recently been observed in cold dark gas in TMC1 and other cores in Taurus (less strongly in HC3N 
there). It was suggested that the predicted fractionation effects had at long last actually been seen, 
albeit in only some very cold gas. We showed here that 13 C fractionation does not occur in two 
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species having a single carbon (CS and HNC) seen toward TMC1 in Taurus; we also observed 
these species toward other positions, finding certainly no clear evidence of fractionation, but 
derivation of an accurate isotopic abundance was frustrated by the difficulty of correcting for finite 
optical depth even in very rare isotopologues. In any case, our observations toward TMC1 make it 
unlikely that any general depletion of 13 C exists in the gas at large outside CO and we discussed 
the implications of this inference in the context of the chemistry of an optically shielded dense gas 
with a normal ratio [0]/[C] > 1 . 

In general, preventing an overall 13 C fractionation requires that the pool of gas-phase carbon 
outside CO resides in the gas mostly as C° rather than C + . This could happen if the [0]/[C] ratio 
is less than unity or if some mechanism is invoked to liberate neutral carbon from CO through 
photodissociation deep inside dark gas but the most interesting possibility is that the TMC1 gas is 
in the so-called high ionization state of a bistable chemical network. In any case, we are left with 
the fact that the overall carbon pool outside CO is apparently not depleted in 13 C even if there is 
no way to prevent 12 C + / 13 C + » 60. For the chemistry this has interesting consequences; reactions 
involving C + see a pool of carbon depleted in 13 C while those with C° do not, and the main pool 
of carbon is in larger amounts of C° that react somewhat more slowly, with a normal 12 C/ 13 C ratio. 

In the future it may be possible to ascertain the overall composition of the carbon pool deep 
inside TMC1 using high spatial resolution observations of sub-mm lines of C I at ALMA, and 
perhaps even to measure the C/ 13 C ratio using the 492.164 GHz line of 13 C that is displaced from 
the main isotope. In the meantime it is important to assess just which chemical species are subject 
to isotopologic abundance anomalies - apparently, a wide variety of tri- and polyatomics hosting 
more than one carbon atom - and under what conditions they arise, and to avoid using strongly 
affected molecles to derive isotopic abundance ratios. 

The Arizona Radio Observatory is operated by Steward Observatory, University of Arizona, 
with partial support through the NSF ATI program. The National Radio Astronomy Observatory 
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Rottach-Egern during the completion of the initial draft of this work and that of Bel Esperance in 
Geneva during the next. We thank the referee for a close reading of the manuscript and a variety 
of interesting suggestions for further discussion. 
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A. Appendix A. The fractionation reaction 



Smith & Adams 



(| 19801) did not give analytic forms for th eir measurem e nts of t he tempera ture 



dependence of the carbon fractionation reaction first cited by 
provided the expressions 



Watson 



et all J 19761) . 



Lisztl (2007) 



k f = 7.64 x 10" y T K " U33 cm 3 s" 1 (T K = 80 - 500K) 



k f = 



1.39 xlO" 9 T K ~ 005 cm 3 s" 1 
1 +exp(-34.8/T K ) 



(T K = 10 - 80K) 



k r = fc/exp(-34.8/T K ) 



for the forward (exothermic) and reverse reaction rates k f and k r 



B. Appendix B. A toy model of 13 C depletion 



We constructed a formal model of the chemistry in which C + is produced by the action of He + 
on a fixed abundance of CO (X(CO)= n(CO)/n(H 2 ) is a parameter) and solved for the densities 
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Fig. 2. — Toy model of the 13 C + , 12 C + chemistry. In each panel the solid line is X(C + ), the (red) 
dashed line is [ 12 C]/[ 13 C]x X( 13 CO)/X( 12 CO) and the (blue) dotted line is X(He + ) (only shown in 
some panels). Standard values are T K = 10 K, n(H 2 ) = 10 4 cnT 3 , ( = 2x 1(T 17 per H 2 , X(CO)= 
8 x 10" 5 . X(M) = 3 x 1(T 7 . X(Oo) = 3 x 10~ 7 and 12 C/ 13 C = 60. The panels of this fisure show the 
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n( 12 C + ), n( 13 C + ) and n(He + ). In addition to the fractionation reaction with CO, C + is destroyed by 
interaction with H 2 , hydrides (CH, NH, OH), 2 , thermal electrons, and charge exchange with a 
single low-ionization metal species denoted "M" and given the atomic properties of silicon. For 



the reactions in Tabl e 4 we used the dipo 



http://www.udfa.net/ (fWoodall et al 



2007). 



e-enhanced rates given on the UFDA06 database website 



The default values of the reactant abundances are given in Table 4 and quoted in Figure 2, 
i.e. n(H 2 ) = 10 4 cm" 3 , T K = 10 K, { H ,= 2x lO^V 1 . For the default electron a bundance 2.4 x 10~ 7 
we used the expression in lOppenheimer & Dalgarnol (|1974|) : alternatively see LMcKeel (|1989|) . The 
default carbon isotopic ratio is Ri2/i3= [ 12 C]/[ 13 C] = 60. The entry for 'M' represents all heavy 
atoms with ionization potentials less than lOev or so and the default is for strong depletion - the 
typical "low metals" case. Nitrogen was assumed to exist in the form of N 2 with a Solar [N]/[C] 
ratio and carbon depletion in the model corresponds to the default value X(CO) = 8 x 10~ 5 . 



For He + 



dn(Ue + )/dt = £ H e"(He) - n(He + )n(H 2 ) V kjXj 



where the fractional abundances of the reactants j are Xj and their reaction rate constants with 
ie + are k,- (li sted in Table 4). For He we take the local galactic disk abundance [He]/[H] = 0.088 



(Balser 



3). 



2006). The cosmic-ray ionization rate of He is £ H e = 1-08 £W2 and ^ H2 = 2 x 10 s .In 



general, direct recombination of atomic ions with ambient electrons is utterly insignificant in this 
context. 



For 12 C + and 13 C + 



dn( u C + )/dt = n( n CO)[n( 13 C + )k f + n(He + )k He . C0 ] 
-n( u C + )[n( 13 CO)k r + n(H 2 ) ^ k._ c +Xj] 



dn( i3 C + )/dt = n( u CO)[n( l2 C + )k r + n(He + )k H e-co] 



,13 
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-n( n C + )[n( n CO)k f + n(H 2 ) ^ k._ c+ Xj] 



In each panel of Fig. 2 the relative abundance X( 12 C + ) is shown by a solid (black) line and 
X( 13 C + ) x R.12/13 (the red dashed line) has been scaled up by the inherent isotopic abundance so 
that the gap between the solid black and dashed red lines shows directly the extent to which 13 C 
is depleted in C + and in the gas outside CO. The effect is generally very large at 10 K, though 
still somewhat below exp(35/T K ). The lower right panel shows the effect of artificially increasing 
the strength of the fractionation reaction, illustrating that depletion of 13 C + might be a factor two 
stronger if competing reactions were less important. 

Although the important chemical effects ahould be incorporated in the default model, its 
13 C depletion (a factor 15) is much stronger than is seen even in CCH (at most a factor 4) and 
substantial 13 C depletion persists at relatively high temperature (the upper left panel). Of course 
none of this is observed. Lowering the C + / 13 C + ratio can be accomplished ad hoc by assuming an 
undepleted metal abundance, by putting all the spare oxygen in a species like O2 or OH that reacts 
rapidly with C + , or by increasing the rate at which C + recombines with H 2 (lower left). Perhaps 



the least ad hoc modification is the case of small X(CO) because C O will eventua 
th e gas phase. This modificatin is basically the effect suggested by 



by 



Liszddl978h . 



ly deplete from 



Watson (1977) and modelled 



C. Appendix C. Is carbon fractionation reversible after formation? 

In the fractionation chemistry the ambient 13 C + is depleted by a factor somewhat smaller 
than exp(35 K/T K ) corresponding to the binding energy difference in CO. In principle, if another 
molecule's chemistry was dominated by a fractionation reaction like that of C + and CO, with an 
energy defect comparable to 35 K and an abundance much less than that of C + (which itself is not 
large) so that the gas contains sufficient amounts of 13 C, that molecule's complement of 13 C could 
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be restored to a degree approaching the abundance ratio inherent in the gas, R12/13. At present 
it is not possible to demonstrate that such in situ fractionation after formation occurs for any 
species beside CO: species either react chemically with C + to form other species or they react too 
frequently with other things beside C + , with the added complication that the rate constants of the 
required fractionation reactions are still unknown even 35 years since the importance of the CO 
fractionation reaction was established. CS is in fact one of the better chemical candidates, with 
an energy difference between CS and 13 CS of 26.3 K, but the fractionation reaction with C + has 
never been measured. 
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Table 1: Species, line frequencies and integrated line temperatures" 



Species 


Frequency 


Cyano Pk 


L1527 


NH 3 Pk 




MHz 


K-km s 1 


K-kms 1 


K-km s 1 


HNC 


90663.56 


1.597(0.020) 


1.850(0.040) 


2.550(0.020) 


HN 13 C 


87090.85* 


0.551(0.023) 


0.263(0.022) 


0.650(0.026) 


H 15 NC 


88865.69 


0.120(0.016) 


0.055(0.015) 


0.086(0.016) 


CS 


97980.95 


0.980(0.011) 


1.495(0.023) 


1.685(0.010) 


c 34 s 


96412.95 


0.362(0.014) 


0.116(0.009) 


0.341(0.005) 


13 cs 


92494.27 


0.121(0.008) 


0.045(0.007) 


0.100(0.008) 


H 2 CS 


103040.28 


0.483(0.012) 


0.102(0.019) 


0.272(0.011) 



H 2 C 34 S 101284.40 0.0218(0.0044) 
H 2 13 CS 99077.84 0.0063(0.0016) 



a Table entries are J T^dv quantities in p arentheses are lcr 



For details of the HN C spectrum see Ivan der Tak et al 



uncer tainty 



(2009) 
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Table 2: Integrated temperature and implied isotopologic ratios 
Species Cyano Pk LI 527 NH 3 Pk 

W(HN 13 C)/W(H 15 NC) 4.59(0.64) 4.78(1.36) 7.56(1.44) 

N(HNC)/N(HN 13 C) fl 54(8)..72(10) 52(15)..69(20) 33(6)..44(8) 

N(HNC)/N(HN 13 C)* 43(6)..57(9) 



W( 13 CS)/W(C 34 S) 0.334(0.024) 0.385(0.066) 0.293(0.028) 

N(CS)/N( 13 CS) C 68(5) 59(10) 77(7) 

N(CS)/N( 13 CS) d 71(5) 

W(H 2 CS)/W(H 2 C 34 S) 22.9(4.2) 

W(H 2 13 CS)/W(H 2 C 34 S) 0.289(0.094) 

N(H 2 CS)/N(H 2 13 CS) e 79(26) 
a for N(HNC)/N(H 15 NC) = 250..330 and r = in the HN 13 C line 

h for N(HNC)/N(H 15 NC) = 250.. 330 and r = 0.6 in the HN 13 C line. 

c for N(CS)/N(C 34 S) = 22.7 and r = in the C 34 S line 

d for N(CS)/N(C 34 S) = 22.7 and t = 0.15 in the C 34 S line 

e for N(H 2 CS)/N(H 2 C 34 S) = 22.7 and r = in the H 2 C 34 S line. 
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Table 3: Other measured and inferred 12 C/ C ratios 



Species 


Cyano Pk 


1521E 


LI 527 


NH 3 


ref 


CCH/ 13 CCH 


> 250* 




> 135 c 




1 


CCH/C 13 CH 


> 170* 




> 80 c 




1 


cc 34 s/ 13 ccs 


230(43) 


>130 






2 


cc 34 s/c 13 cs 


54(2) 


51(4) 






2 


c 3 34 s/ 13 cccs 


> 191 


> 107 






2 


HC 3 N/H 13 CCCN 


79(11) 








3 


HC 3 N/HC 13 CCN 


75(10) 








3 


HC 3 N/HCC 13 CN 


55(7) 






45(6) 


3 



a Lower limits are 3cr, rms in parentheses are lcr 
b C 13 CH/ 13 CCH = 1.6 ±0.4(3(r) 
c C 13 CH/ 13 CCH = 1.6 ±0.1(3(r) 



R eferences: 1) 
2) 



3) 



Sakai etal 



SakaietaL|(2010) 



(2007) 



Takano et all (1 1998) 
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Table 4: reactants and their rate constants with He + and 



Reactant Relative abundance 



w/He + 



3 -I 

cm s 



w/C 4 



s -1 
cm s 



H 2 
CO 

e 

M 

CH 

OH 

2 
H 2 

^2— 



1 

8.0 x 10" 5 
2.4 x 10~ 7 
3.4 x 10~ 7 
1.0 x 10" 7 
3.0 x 10~ 7 
1.0 x 10~ 7 
l.Ox 10~ 7 

-5 



7.2 x 10" 15 +3.7 x 10" 14 exp (-35/T K ) 4.0 x 10" 16 (300/T K ) 



\0.20 



1.6 x 10 9 

3.3 x 10 9 
1.6 x 10' 9 V(300/T K ) 
1.0 x 10' 9 V(300/T K ) 

l.Ox 10~ 9 



4.8 x 10- 10 V(300/T K ) 
1.6 x 10 9 



see Al 
7.2 x 10~ 12 (300/T K )°- 83fl 
1.5 x 10~ 9 
3.8 x 10- 10 V(300/T K ) 
7.7 x 10- 10 V(300/T K ) 
1.0 x 10' 9 



2.7 x 10 9 V(300/T K ) 



a Wolfire et al. (2008), includes dielectronic recombination 



